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ABSTRACT 

We develop a model for stochastic pre-enrichment and self-enrichment in globular clusters (GCs) 
during their formation process. GCs beginning their formation have an initial metallicity determined 
by the pre-enrichment of their surrounding protocloud, but can also undergo internal self-enrichment 
during formation. Stochastic variations in metallicity arise because of the finite numbers of supernova. 
We construct an analytic formulation of the combined effects of pre-enrichment and self-enrichment 
and use Monte Carlo models to verify that the model accurately encapsulates the mean metallicity 
and metallicity spread among real GCs. The predicted metallicity spread due to self-enrichment alone, 
a robust prediction of the model, is much smaller than the observed spread among real GCs. This 
result rules out self-enrichment as a significant contributor to the metal content in most GCs, leaving 
pre-enrichment as the viable alternative. Self-enrichment can, however, be important for clusters with 
masses well above 10 6 M , which are massive enough to hold in a significant fraction of their SN ejecta 
even without any external pressure confinement. This transition point corresponds well to the mass at 
which a mass-metallicity relationship ( "blue tilt" ) appears in the metal-poor cluster sequence in many 
large galaxies. We therefore suggest that self-enrichment is the primary driver for the mass-metallicity 
relation. Other predictions from our model are that the cluster-to-cluster metallicity spread decreases 
amongst the highest mass clusters; and that the red GC sequence should also display a more modest 
mass-metallicity trend if it can be traced to similarly high mass. 

Subject headings: globular clusters: general — galaxies: abundances — galaxies: star clusters - 
galaxies: formation — galaxies: evolution — methods: analytical 



1. INTRODUCTION 

Globular clusters (GCs) are believed to form in sin- 
gle bursts of star formation within protocluster clouds, 
producing bound clusters of 10 4 - 10 7 stars which (in 
most cases) share a single age, metallicity, and element 
abundance pattern. A near-universal feature of the GCs 
in any one large galaxy is that they typically follow a 
bimodal metallicity distribution, a conclusion now based 
solidl y on large statistical samples from m a ny galaxies 
/e.g. lLarsen et all 120011 lHarris et al] 120061 : iPeng et all 
120061 120081 ). The more metal-poor (blue) clusters av- 
erage [Fe/H] ~ -1.5 with a dispersion cr[Fe/H] ~ 0.25, 
and metal-rich er (red) ones avera ge [Fc/H] ~ —0.3 with 
cr[Fe/H] ~ 0.4 (lHarris et al.ll200ll . 

A more recently discovered second-order feature in the 
metallicity distributions, that seems particularly to af- 
fect the metal-poor component, is a mass-metallicity re- 
lation (MMR) : the blue clusters become slightly but sig- 
nificantly mo r e met a l-rich at progressive ly higher mass 
(lHarris et alJ 120061 : IStrader et at] 120061 : IMieske et ail 
2006). No such correlation is seen along the metal-richer, 
red GC sequence. This feature, also colloquially called 
a "blue tilt" , is as yet poorly understood, and the vari- 
ous empirical characterizations of the effect are not yet 
in mutual agreement, even including the extreme claim 
that it is only an artifact of measurement (|Kundul [2008) . 
Some discussions of the MMR in a few individual galaxies 
claim to find a systematic correlation of cluster metallic- 
ity (color) with mass (luminosity) that spans the entire 
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GC luminosity range (IStrader et all 120061: ISpitler et alJ 
l200aiWehner et alJl2008l: iSpitler et al.ll2008l) . However, 
the observational studies based on the biggest statistical 
samples (thousands of GCs collected from many large E 
galaxies) , and carried out with the best photometric mea- 
surement techniques (with careful attention to aperture- 
size corrections determined from convolutions of the GC 
profiles with the point-spread function of the image) find 
that the slope of the MMR becomes most noticeable 
for the top end of the GC mass range, M j> 10 6 Mq 
(jHarris et al.ll2006l : IMieske et al.ll2006l : (Harrisll2009l ) . For 
the vast majority of lower-mass clusters, their mean color 
stays virtually constant and it is debatable whether or 
not any significant MMR trend exists for M < 10 6 M Q . 

The fact that the MMR becomes prominent only at the 
high- mass, high-luminosity end of the GC sequence easily 
explains why it was not discovered until recently. Mas- 
sive GCs are rare, and for galaxies with relatively small 
GC populations (such as the Milky Way, any dwarfs, 
or any non-giant ellipticals), the GC distribution simply 
cannot be traced to high enough levels for any systematic 
change in cluster color with luminosity to be noticed. In 
the Milky Way, for example, only u) Centauri is luminous 
and massive enoug h to be clearly in the "MMR regime" 
(|Harris et al.l [2006) . To make the picture additionally 
puzzling, one or two giant ellipticals seem to have no 
discernable MM R at any luminosity level (most notably 
NGC 4472; see IStrader etTall l2?M IMieske et al.l 127)061 ) 
and it is not out of the question that the amplitude of 
the effect may differ from one galaxy to another. 

Two other lines of observational evidence have recently 
emerged to hint more strongly that the systematic prop- 
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erties of GCs start to undergo important changes at 
around the million- Solar-mass point. One of these is 
their characteristic radii (usually described by their half- 
light or effective radius fh, which is relatively immune 
to dynamical evolution). Empirically, r^ stays relatively 
uniform at r^ ~ 3 pc for luminosities L < 10 6 Lq (cor- 
responding roughly to M < 2 x 10 6 M Q ), but follows a 
steeper relation rou ghly approximated b y r h ~ L 1 / 2 at 
higher luminosities (iBarmbv et ail 120071 : iReikuba et al.l 
120071: lEvstigneeva et al.ll2008t) . This changeover in their 
structural-parameter plane may smoothly link the most 
massive GCs with other objects such as UCDs (ultra- 
compact dwarfs) and the compact nuclei in dwarf ellip- 
ticals. 

Lastly, evidence gathered from detailed color- 
magnitude studies and abundance measurements of 
GCs within the Milky Way (e.g . iBedin et all 12004 
Piotto et all 120071: IWlanova et all l2007t iMilone et al l 
20081: Uohnson et all l200l is now showing that the 
most massive GCs (again, those typically above about 
10 6 Mq) are often not the clean "single stellar popu- 
lation" that is classically associated with GCs; instead, 
they can show two or more distinct sequences in their 
color-magnitude diagram that are indicators of more 
than one generation of stars with different chemical abun- 
dances. 

GC formation models to explain and link all these in- 
triguing characteristics of GCs at different masses are 
still in rudimentary stages. In this paper, we develop 
a new quantitative model to address particularly the 
metallicity distributions of GCs and their systematic 
change with mass. Our model assumes that a globular 
cluster in the process of formation can be pre-enriched 
(that is, its protocluster gas will have an initial heavy- 
element abundance produced by any generations of stars 
that happened previously), and also self-enriched (that 
is, its own first rounds of supernova can further increase 
its mean metallicity). Our model is intended to explore 
which of these factors will be dominant in a given situa- 
tion, and whether the combination is capable of describ- 
ing the metallicity distribution for real GCs. 

Because all timescales associated with the formation 
and evolution of stars decrease strongly with increasing 
stellar mass, it is possible for the first high-mass stars 
in a protocluster to form, live their entire lives, and ex- 
plode as supernovae (SNe) while the lower-mass stars 
are still forming. The metals produced from high-mass 
stars within a protocluster may therefore be incorporated 
into the lower-mass stars in the same starburst (self- 
enrichment). In this case, the mean metallicity of the 
final cluster will reflect several factors including the ini- 
tial mass function (IMF) of star formation, the efficiency 
with which the protocloud gas is turned into stars, and 
the ability of protoclusters to hold on to the metals pro- 
duced by its SNe (which will be primarily a function 
of the cloud mass). Because the number of SNe in a 
protocluster may be subject to small-number statistics, 
the amount o f metal enrichment can be quite stochastic 
(|Cavrellll986f) . resulting in a metallicity spread that itself 
will be a f unction of cluster mass. 

Recentlv lStrader fe Sm ith (2008]) have published a for- 
mation model for GCs also based on the idea of self- 
enrichment. We compare their discussion with ours in 



more detail in a later section below. We believe that 
key differences in the assumptions regarding the proto- 
cloud structure and formation timescales make our model 
more physically realistic, and as will be seen below, make 
it possible to identify the intriguing "transition point" 
around 10 6 Mq in a more natural way. 

2. STAR FORMATION AND EVOLUTION TIMESCALES 

Previous descriptions of GC and star-cluster formation 
usually make the zero-order approximation of an "in- 
stantaneous starburst" in which the entire stellar IMF 
springs into place at once. At some level this must be an 
oversimplification: individual stars have formation times 
that depend on their mass, and the starburst epoch as 
a whole will be stretched out over a comparable length 
of time. Our model of self-enrichment is based instead 
on the picture that many or most of the high-mass stars 
form, live, die, and eject their newly formed heavy el- 
ements into the surroundings before the low-mass stars 
have finished assembling a significant fraction of their 
mass. This scenario, too, must at some level be an ide- 
alization, but there is now significant evidence on both 
the observational and theoretical sides to justify it as a 
useful first-order approximation. 

In the context of globular cluster self-enrichment, 
"high mass" stars are those that contribute significantly 
to the total metal production from the stellar popula- 
tion while "low mass" stars are those that are still vis- 
ible today in a Milky Way GC. The majority of metals 
are produced in supernovae with progenitor masses of 
at least 20 Mq (see Figure [lj, while the upper mass 
limit for "low mass" stars relevant here will be the main- 
sequence turnoff mass of typical GCs in large galax- 
ies, Mto IMq. Therefore, the important stellar 
timescales are the combined formation and evolution 
timescale of stars of at least 20 Mq compared to the 
formation timescales of stars of less than 1 Mq . 

The timescales associat ed with high - mass star forma- 
tion are extremely short. IChurchwelll (|2002t ) has shown 
that the dynamics of molecular outflows in massive star 
forming regions require very high accretion rates that 
imply formation timesc ales of 10 4 - 10 5 yr (see also 
iMcKee fc Ostri kcr 2007, who infer that the formation 
times for massive stars are < 1 Myr). Once formed, 
the ev olutiona r y lifet ime of high mass stars is similarly 
short. iHirschil (|2007f ) calculates lifetimes of 10 Myr for 
20 Mq star s down to < 6 Myr for s tars of 40 Mq or 
greater, and lRecchi fc Danzigeri (|2005l ) model combined 
effect of metal production and stellar ages and find that 
most metals are released within the first 5 Myr after the 
first 100 Mq star dies (see their figure 6). The timescale 
for the supernova ejecta to cool is similar, ~ 7 Myr 
(jMcKee fc Ostrikerlll977h . Therefore, even though stars 
with masses as low as 8 M Q may explode as supernovae, 
the key timescale for enrichment is much shorter than 
the ~ 50 Myr lifetime of those stars. In short, a reason- 
able estimate appears to be that the high-mass part of 
the IMF can 'seed' its surroundings beginning around 5 
Myr after formation. 

On the other hand, the formation times for low mass 
stars are much longer. The formation timescale should 
be approximately the Kelvin-Helmholz time, which rises 
to t K h > 10 7 yr for M < 1 Mq (see figure 4 of 
IZinnecker fc Yorke! I2007T) . Therefore, even if all stars 
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in a given starburst begin to form simultaneously, the 
heavy elements from the high mass stars will be ejected 
into the protoclustcr cloud while the low mass stars are 
still collapsing. In fact, simulations that couple hydro- 
dynamics with radiative transfer suggest that feedback 
from high mass stars critical ly shapes the low mass IMF 
itself (|Krumholz et al.|[2007t ), providing direct evidence 
that feedback from high mass stars is important early in 
the formation of lower mass stars. 

These comments still do not include the internal 
structure of the starburst as a whole. The turbulent 
medium in which star formation occurs is theoretically 
expected to result in clumpy star formation, as is seen 
observationally in dense star-forming regions toda y (e.g. 
iFeigelson fc Townslevi[200l iTownslev et al.ll2006l) . both 
in the spatial structure of the protocluster and in its tem- 
poral evolution over the star formation episode. These 
studies indicate that the total durations of starbursts in 
massive, dense protocluster clouds such as the ones we 
are particularly interested in here are plausibly in the 
range of 10 — 20 Myr. Intervals this long give the low- 
mass stars that are formed throughout the burst addi- 
tional time to incorporate metals ejected by the higher- 
mass stars formed at the onset of the burst. While we 
will not address the effects of this internal structure on 
the process of cluster self-enrichment any further, it may 
be expected to introduce star-to-star metallicity disper- 
sion within the cluster; inhomogeneities of this sort have 
indeed been detected in some massive GCs (see the ref- 
erences cited above). 

3. STOCHASTIC MODEL OF SELF-ENRICHMENT 

We begin by analyzing the simplest case, where 
the metal content of clusters comes entirely from self- 
enrichment, in order to address whether it can be the 
dominant contributor. We expand this analysis to in- 
clude contributions from both self-enrichment and pre- 
enrichment in § [5] 

3.1. Mean metallicity 

We assume that a globular cluster forms from a molec- 
ular cloud of mass M c , with a star formation efficiency of 
/*. The stars are distributed according to an IMF £(m) 
of the form 



_ dn 
C(m) = d^ 



Am c 



(1) 



between a minimum mass M m i n and a maximum mass 
M msx (assumed to be 100 M ). We assume a Salpeter 
slope of a — —2.35, which is an excellent fi t to the ob- 
served IMF at m > 1 M Q (|Chabrierll2003] ). Although 
the low mass IMF has large uncertainties, it is known to 
turn over and therefore contain only a small fraction of 
the total stellar mass. Our derivation only requires that 
(a) the high-mass slope of the IMF and (b) the total 
fraction of mass contained in high mass stars are cor- 
rect. We may therefore assume a Salpeter IMF with an 
appropr iate low mass c utoff M min without loss of gen- 
erality. IChabrieil (|2003f ) finds that the mass fraction of 
stars with m > 9 M Q is 0.20; for our adopted IMF, this 
implies M m i„ = 0.30 M Q . 

The normalization constant A can be found from the 
condition that the integrated mass over the IMF is equal 




40 60 80 
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Fig. 1. — Cumulative fraction of metals produced by supernovae 
with progenitor star masses less than or equal to the given mass. 
Almost 80% of the metals come from stars with initial masses of 
at least 20 M Q . 

to the total stellar mass Mqc (= f*M c ): 



A = 



f*M c (a + 2) 



a+2 



M n 



a+2 ' 



(2) 



We assume that all stars above a critical threshold 
Msn ~ 8 Mq explode as SNe. The mass of metals 
mz ejected by core collapse supernovae from stars of 
various progenitor masse s m have been ca l culate d by 
IWooslev fe Weaverl (j!995h and lNomoto et alj (| 19971) . As- 
sessment of their data shows that they are reasonably 
approximated by the relation 



mziyri) = (B + Cm)m, 



(3) 



where B = 1.18%, C = 0.548%, and all masses are in 
units of M Q . Thus, for example, a 20M Q star would 
return to its surroundings about 2.4M of enriched gas 
in all heavy elements combined. We assume that these 
yields do not depend on the metallicity of the progenitor 
stars, an assumption we justify in § 16.1.21 The total 
amount of metals produced by SNe within the cloud is 
therefore 



M z = 



mz(m)^(m)dm 



(4) 
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a + 3 
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a+3 



M SN 



SN 



a+3\ 



a+2 



) + 



(5) 



The relative importance of stars of different mass can be 
seen in Figure [1] where we plot the cumulative fraction 
of total metals contributed by supernovae as a function 
of the progenitor star mass (i.e. the plotted value is 

f( M ) = JM SN m z(m)£(m)dm/ m z {m)£(m)dm). 
Although the IMF rises to low masses, the yield per star 
drops almost proportionally and so stars with masses 
near Mqn do not contribute appreciably to the total 
metal production; almost 80% of the metals come from 
stars with initial masses of at least 20 M Q . 

We assume that a fraction fz of the metals produced 
stays within the the cloud and ends up incorporated in 
the formation of the lower-mass stars. For the first stage 
of our discussion, we assume that fz is a constant free 
parameter; we will deal with the more realistic case where 
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fz depends on the depth of the potential well in § [5j The 
total metallicity Z c of the gas from which the stars are 
formed is then 

„ fzM z , . 

Zc = (6) 

with Mz taken from equation (J5]). The resulting metal- 
licity of the low-mass stars that we see in the cluster 
today, many Gy later, is then 



logZ c /Z = 0.38 + log/, f z , 



(7) 



where we adopt Zp, = 0.016 JG revessc k SauvallH99l 
IVandenBerg et al1l2007f) although the exact value is not 
critical to our argument. We note that our parameter fz 
is essentially similar to the "effective yield" y e ff in sim- 
ple models of chemical evolutio n (e.g. iPagel k PatchettJ 
[19751: iBinnev k MerrifjeldlflOOSh . 

There is no reason to expect either the star forma- 
tion efficiency /* or the metal retention efficiency fz 
to b e near unity. Observations indicat e that /* ~ 
0.3 |Lada et all 11984 lLada k Ladall2003t iMarks et all 
2008). Though fz is less well constrained by direct ob- 
servation, we can use the above derivation to constrain 
the product f*fz- Because the amount of self-enrichment 
sets only a lower bound to the total GC metallicity, 
log Z c / Zq should not be greater than the mean metal- 
poor GC metallicity of [m/H] ~ —1.25, which implies 
log/*/z < -1.63 and thus (if /* - 0.3) f z < 0.08. 

Furthermore, /» and fz may in principle vary as a func- 
tion of cloud mass. For example, the fraction of metals 
that remain bound to the cloud should increase with the 
depth of the gravitational potential well, and therefore fz 
may increase with mass. This by itself would introduce 
a positive mass-metallicity relation in the same sense as 
the observati ons listed above. We d iscuss this in more 
detail in § |5J IStrader k Smith! (|2008f) make an argument 
that is basically similar. 

3.2. Metallicity spread 

Next, we extend the above derivation to determine the 
expected cluster-to-cluster dispersion in the GC metal- 
licity distri bution, due t o the stochastic nature of self- 
enrichment (jCayrel 1986|). This will turn out to yield an 
interesting constraint on the relative importance of pre- 
and self-enrichment. 

First we note that the mean number of supernovae in 
a given GC, N~sn, is: 



Ni 



SN : 



Mi. 

M S N 

A 



1 



£{m)dm 
(M, 



a+l 



Mi 



SN 



a+l\ 



(8) 
(9) 



The variance in the mass of metals expelled by a single 
SN is: 

(10) 



(m|) - {m z } 



where 
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and 



(mz) 



Ml 
^sY 



(12) 



(13) 



with Mz taken from_ equation ((5|). 

Adding together TVsn supernovae, the variance in the 
total metal mass Mz is 



(14) 

and the relative scatter in the cloud metallicity Z c is: 

&z _ °M Z 
~ ~M~z' 

The original cloud mass M c is not observable, but the 
mass in stars Mqc is. Combining equations @, (O, and 
© - (HU), we find: 



(15) 



^ = 0.059 ( 

z c Vio 5 m q 



-1/2 



(16) 



This formulation reveals a remarkable result: the pre- 
dicted width of the metallicity distribution is a function 
of only the stellar mass in the globular cluster, indepen- 
dent of the free parameters /* and fz- It is a direct 
consequence of the fact that higher cluster mass yields 
larger numbers of SNe and so reduces the statistical fluc- 
tuations in the enrichment, regardless of how much gas is 
lost. This relationship, including normalization, is there- 
fore a robust testable prediction of the model. 

3.3. Monte Carlo realization 

We next generate a Monte Carlo realization of this pro- 
cess to demonstrate the stochastic effects more explicitly. 
We begin by randomly sampling cloud masses M c from 
a mass function for protocluster molecular clouds of the 
form 

dN 



dM r 



oc M: 



(17) 



between a minimum mass Mq and M — > oo 
(jHarris k Pudritzl |1994| estimate a slightly shallower 
slope for the mass function, ~ —1.6, but our results are 
entirely insensitive to its value). 

For each cloud, we calculate the mean number of stars 
expected, N: 



N = 



t;(m)dm 



UM c {a + 2) M n 



a+l 



Mrr 



a + l 



(a + 1) M B 



a+2 



Mr, 



a+2 ■ 



(18) 



(19) 



and round it to the nearest integer. We then randomly 
sample N stars from the IMF and find the stellar mass of 
the cluster, Mqq, by taking their sum. Because the input 
cloud mass M c and the mass obtained from a random 
sampling of the IMF may not agree, we reassign the cloud 
mass to M c = Mqc //* 3 - For each star with mass greater 

3 We have confirmed that the distribution of reassigned cloud 
masses also follows equation H71 . 
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o) f.= 1 f z =1 
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f. = 0.1 f z =1 

». 
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c) 


f.= 1 f z = 0.1 - 

f. 


1.0 -0.5 0.0 0.5 
log Z/Z s 



-1.0 -0.5 0.0 0.5 
log Z/Z 

Fig. 2. — The stellar mass of globular clusters, Mqq is plotted 
against their metallicity, logZ c /ZQ. Panels (a), (b), and (c) corre- 
spond to the models /* = 1, fz = 1; /* = 0.1, fz = 1; and /* = 1, 
fz = 0.1 respectively. The mean metallicity from equation (0 
is shown as t he red solid line, and the predicted dispersion from 
equation ||16|I is shown by the blue dashed lines. Dots indicate the 
Monte Carlo points. 
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Fig. 3. — Number of supernovae per cluster as a function of clus- 
ter mass, for the Monte Carlo realization (dots) and the expected 
mean from the equation J9] red line). The expected Poisson dis- 
persion is shown as the blue dashed lines. 

than Msn, we calculate the mass of metals it produces 
using equation ([3|), and finally calculate the metallicity 
of the cloud by summing up the contributions from each 
SN and multiplying by fz/M c . 

The results of 1500 Monte Carlo models with a mini- 
mum cloud mass of M = 10 4 //* Mq (so that the mini- 
mum Mqq is constant), along with the results from the 
analytic derivations, are given in Figure[2] It can be seen 
that the analytic results and the Monte Carlo samples 
are completely consistent, and that the free parameters 
/* and fz only affect the mean metallicity and not the 
scatter at a given globular cluster mass. 

We plot the number of supernovae per cluster in Fig- 
ure [31 The average is ~ 1 SN for every 100 Mq of stars 
formed. Thus, for example, a 10 5 Mq globular cluster - 
that is, an object that is thought of as a classically "nor- 
mal" GC - will experience on average 1000 supernovae 
during its formation. Note that the spread in cluster 
metallicities is not simply the Poisson error of the num- 
ber of supernovae, but is also due to the spread in pos- 
sible metal contributions of each individual supernova. 
For example, the relative Poisson error on 1000 super- 
novae is ~ 0.03, while the actual metallicity spread for 
clusters of Mqc = 10 5 Mq is a factor of two larger. 

Lastly, we emphasize that the calculations so far as- 



10" r 



— 10° 



10 4 : 




-3 -2-10 1 

log Z/Z 

Fig. 4. — The red solid line and blue dashed lines show the mean 
metallicity and its spread predicted from our pure self-enrichment 
model for /* = 0.3 and fz = 0.08, while the black dots represent 
1500 Monte Carlo realiz ations. The gr een diamonds represent the 
observational data from Harris (1996) for Galactic globular clus- 
ters, and have been plotted assuming log Z/Zq = [Fe/H] + 0.25. 

sume a constant gas retention fraction fz, so that the 
mean self-enrichment is independent of cluster mass. 
This assumption is rather obviously unrealistic, since at 
small enough protocluster masses the first few super- 
novae will remove all the gas and fail to generate any 
enrichment. We deal with this feature in Section 4 be- 
low. 

4. COMPARISON OF SELF-ENRICHMENT MODEL WITH 
OBSERVATIONS 

In Figure [H we compare the metallicities and stellar 
masses predicted by our self-enrich ment model to the ob- 
served Milky Way globular clusters (|Harrislll996h . When 
comparing to observations, we assume 



log Z/Z G = [Fe/H] + 0.25, 



(20) 



with a scatter of 0.1 dex (e.g. [G ratton & Qr tolanil 
1989tlShetrone et aU200lHPritzl et al. 2005; Kirb v et all 



20081) . We adopt /* = 0.3 and f z = 0.08 in order that 
the predicted mean metallicity matches the mean metal- 
licity of the observed metal-poor clusters. 

The observed metal-poor clusters have a metallicity 
spread of cr[Fe/H] sw 0.27 with no obvious dependence 
on cluster luminosity. This corresponds to a relative 
metallicity spread of <7 Z /Z = cr[Fe/H] In 10 = 0.62. 
In contrast, the self-enrichment model predicts a rela- 
tive metallicity spread of only ~ 0.02 for clusters with 
Mgc — 10 6 Mq, rising to ~ 0.2 for very low-mass clus- 
ters with Mqq — 10 4 Mq. In other words, the observed 
cluster-to-cluster metallicity spread is over an order of 
magnitude larger than the spread predicted by our self- 
enrichment model at most cluster masses. 

If self enrichment were the dominant source of metals 
in the metal-poor globular cluster population, then the 
observed and predicted scatter should be similar, and 
should be correlated with cluster luminosity. The dra- 
matic discrepancy between the observed and predicted 
spread, along with the absence of any observed rela- 
tion between cluster luminosity and metallicity spread, 
strongly suggests that self- enrichment is not the domi- 
nant source of the heavy elements in most metal-poor 
globular clusters. We conclude that these clusters must 
have taken their heavy-element compositions from the 
protocluster clouds before star formation began ("pre 
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enrichment"), and that the observed scatter reflects the 
variety of environments in which clusters formed. 

5. THE METAL RETENTION EFFICIENCY AND A 
PREDICTED MMR 

In § 13.11 we derived an upper limit of fz < 0.08 for the 
metal retention efficiency based on the mean metallic- 
ity of metal-poor globular clusters. Moreover, the upper 
boundary fz ~ 0.08 corresponds to clusters whose met- 
als are entirely contributed by self-enrichment, while we 
concluded in §|4]that self-enrichment cannot be the dom- 
inant source of metals. It is natural to ask whether such 
low values for fz are physically plausible. In this sec- 
tion, we develop a simple energetics argument to obtain 
a rough estimate for fz and explore the consequences of 
the mass- dependent metal ret e ntion that naturally arises 
(see also IDekel fc Silld fl98l IDekel fc Wool 120031 who 
use a similar energetics argument to examine supernova- 
driven gas removal from dark matter halo-enshrouded 
dwarf galaxies). 

For simplicity, we assume that the cloud can be ap- 
proximated by a truncated singular isothermal sphere 
(TSIS), with a density distribution of the form: 




p(r) 



where r t is the truncation radius, 
the cloud is given by 



(21) 



if r < r t , 
if r > r t , 
The potential within 



$(r) 



GM C 

n 



In r. 



(22) 



The total kinetic energy imparted into the cluster gas 
is the sum of the energy injected by each supernova. In 
principle, the important quantity is the number of super- 
novae that explode within a cooling time, and therefore 
the supernova rate, rather than the total number, is im- 
portant. However, in practice the timescale for high-mass 
star formation is expected to be < 1 Myr (see discussion 
in § [2]) , sufficiently shorter tha n the ~ 7 Myr in which a 
supernova bubble remains hot (|McKee fc Ostrikeriflimf) 
that all of the cluster supernovae can be approximated 
as simultaneous. 

We next assume that all of the supernova energy is 
thcrmalizcd within the cluster. If each supernova imparts 
an energy Esn, then the specific kinetic energy K of the 
gas after thermalization is 



K 



NsnEsn 



(23) 



If the kinetic plus potential energy of the gas is greater 
than the potential energy at the edge of the cloud, then 
the gas can escape. This occurs for gas beyond an "es- 
cape radius" , r osc : 



$(r csc ) + X = $(r t ). (24) 
.-.r csc = r te -*s»E SN r t /GMl (25) 

Because the enclosed mass in a TSIS is proportional to 
radius, the fraction of gas within r csc is simply r csc /r t . 
Assuming that the metals are fully mixed, the metal re- 
tention fraction is equal to the gas retention fraction, and 
is equal to 

f z — e -N 3N E 3N r t /GM^ ^6) 
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Fig. 5. — Predicted metal retention efficiency fz as a function of 
protocluster cloud mass (bottom axis) and final cluster mass (top 
axis, assuming a star formation efficiency of /„ = 0.3). The solid 
line is for the fiducial TSIS case, while the other line styles arc the 
predictions for different protocluster density profiles (see § I6.4I I. 

TVsn scales with M c , and is well described by 



A^SN 



M, 



f*Mc 

10 2 M n 10 2 M G 



GC 







Equation (|26[) can therefore be recast as 



fz 



exp 



10 2 MnGM c 



(27) 



(28) 



In Figure [5l we have plotted this relation assuming 
£ SN = 10 51 erg, /* = 0.3 and r t = 1 pc. Note 
that our adopted radius is smaller than the common 
present-day cluster radius of 3 pc due to the expan- 
sion of protoclusters upon remova l of their gaseous enve- 
lope (|Baumgardt fc Kroupall200"7t ); if we adopted a larger 
value for r t , the curve in Figure [5] would maintain its 
shape but move to proportionally to the right, i.e. to 
higher mass. It is apparent that fz -C 0.08 for most 
cluster masses, as required by our argument that self- 
enrichment is unimportant. 

Another way to express this is in terms of the char- 
acteristic cloud mass at which metals become efficiently 
retained, M™ tain , which we define to be when the argu- 
ment of the exponential reaches —1: 



M. 



retain 



E S N.f*r t 
10 2 M G G' 



(29) 



For the values given, M™ tam 
the observable stellar mass, 



4 x 10 7 M H . In terms of 



M rctai„ = J^retain 



EsNf^t 



(30) 



and has a value of Mgp™ 



10 2 M G" 

1 x 10 7 M©. This is at 
the upper range of GC masses, and signals that self- 
enrichment may begin to become important for the most 
massive clusters. 4 

The assumption of complete mixing of the enriched gas 
likely results in an overestimate of fz (or equivalcntly, 

4 We note that the mass Mqq as we define it here is the mass 
just after star formation finishes up. The observed mass of a GC 
as we see it today is only a lower limit to Mqq because of con- 
tinual mass loss from SNe, stellar winds, and ongoing dynamical 
evaporation of stars and tidal trimming. For the massive clus- 
ters that we are particularly interested in, t his mass loss over sev- 
eral Gy may be as much as a factor of two (Vcspcrini ct afl l2003l : 
Baumgar dt fc Kroup a 2007; Wcidn cTet al.11200711 . 
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Fig. 6. — Predicted globular cluster mass-metallicity relation- 
ship assuming a constant pre-enrichment of log Z/Zq = — 1.4 and 
self enric hme nt according to our model, with fz derived from 
equation H28I I. The dotted, dashed, and solid lines indicate the 
metallicity due to pre-enrichment, self-enrichment, and their sum 
respectively. The shaded region indicates the predicted scatter, 
assuming that the relative scatter in pre-enriched metallicity is 
<r[m/H] = 0.27. We have adopted /„ = 0.3; increasing (decreas- 
ing) /* will increase (decrease) the threshold mass at which self- 
enrichment kicks in by a factor of /J while simultaneously increas- 
ing (decreasing) the saturation metallicity at very large Mqq by a 
factor of /*. 

an underestimate of M retam ) compared to the real case 
where the cluster wind contains a disproportionate frac- 
tion of the metals produced in the supernovae. The as- 
sumption of complete energy thermalization may either 
result in an overestimate of fz because it dilutes the su- 
pernova energy by the full cluster mass rather than by 
a smaller fraction of the material, or an underestimate 
of fz because it assumes that no energy is lost to ra- 
diative c ooling. Using more deta iled dynamical consid- 
erations, iParmentier et all (|1999f ) argue that protoglob- 
ular clouds of significantly lower mass can contain their 
supernova ejecta than would be expected from simple 
energy considerations (because of their two-phase model 
where the enriched g as comes from a prev ious genera- 
tion of SNe; see also iMorgan &; Lakel 119891 ), indicating 
that our derivation may underestimate fz (overestimate 
M retain ). While the relative magnitudes of these effects 
are difficult to assess without full blown hydrodynamic 
simulations, their opposite signs reassure us that our re- 
sults are unlikely to be incorrect by orders of magnitude. 
Given the approximations inherent in our derivation, we 
caution against over-interpretation of the detailed form 
of fz, but expect that the general features — particu- 
larly a transition from no self-enrichment at low masses 
to significant self-enrichment above a threshold GC mass 
around Mqq* u — are robust. 

In Figure[6l we have plotted the mass-metallicity rela- 
tion that would be expected if metal-poor GCs are pre- 
enriched to ([Fe/H]) = -1.65 (\ogZ/Z Q = -1.4), and 
then further self-enriched according to our model, noting 
again that our adopted supernova yields are independent 
of the pre-enrichment level. The dotted line indicates 
the pre-enriched metallicity, the dashed line indicates the 
metals added by self enrichment, and the solid line is the 
total predicted metallicity. For these sets of parameters, 
we see that self-enrichment is predicted to begin con- 
tributing at a stellar mass of ~ 4 x 10 6 M Q , very nearly 
where the observed mass-metallicity relation begins. The 
shaded region indicates the predicted scatter, where we 




1.0 1.5 2.0 2.5 

(B-0o 



Fig. 7. — The lines and shaded regions denote the predictions of 
our combined pre-enrichment plus self-enrichment model for pre- 
enrichment values of [m/H] = —1.6 with scatter cr[m/H] = 0.27 
(left sequence) and [m/H] = —0.4 with scatter <r[m/H] = 0.40 
(right sequence), translated into the observational color-luminosity 
space. Line styles are as in Figure \E\ We have ove rplotted the 
observations of giant elliptical GCs from Harris (2009) for compar- 
ison. 

have assumed that the relative scatter in the pre-enriched 
metallicities is cr[m/H] = 0.27, as required to fit the ob- 
served spread among Milky Way GCs. Because the rela- 
tive scatter due to stochastic self-enrichment is very small 
at high masses, the total scatter is entirely dominated 
by the scatter in the pre-enrichment values. Therefore, 
the metallicity scatter in absolute units is constant, and 
the relative scatter becomes smaller at higher mass. Al- 
though the details of the mass-metallicity relation are 
sensitive to the rough estimate of fz, and therefore 
should not be over-interpreted, it is a clear prediction 
of our model that the metallicity scatter in dex should 
decrease as one goes up the mass-metallicity relation, 
preserving an approximately constant absolute scatter 
o-z- 

We directly compare the observed and predicted 
MMRs in Figure [71 usi ng the data fo r 6 giant ellipticals 
recently reme asured by | Harrisl (1 2009) from the original 
data sample of lHarris et al.l (|2006h . Here, the metal-poor 
and metal-rich sequences assume mean pre-enrichments 
to [m/H] = —1.6 and [m/H] = —0.4 with scatters of 
cr[m/H] = 0.27 and <j[m/H] = 0.40 respectively. These 
mean metallicities are chosen deliberately to match the 
blue and red GC sequences. We have converted the pre- 
dicted metallicities and stellar masses into observed ab- 
solute magnitudes and colors assuming 

Mi = 3.94 -2.5 log— , (31) 

Lq 

(B-I)o = 2.158 + 0.375 log — (32) 

Z© 

(]Harris et al.| [2006). We adopt an /-band mass-to-light 
ratio of M/L i = 1.5 M^/Lp), co rresponding to M/Ly — 
2.0 Mq/Lq (jMcLaugh lin 2000) and a mean GC color 
V-I= 1.0 

The model for the metal-poor clusters shows the same 
general features as the observations - a constant metal- 
licity and scatter at most luminosities, with an MMR for 
the most luminous clusters. The observed MMR begins 
at a lower luminosity than the predicted MMR; this may 
be related to our overestimat e of M retain compared t o 
the more detailed dynamics in IParmentier et "all £l999), 
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Fig. 8. — As in Figure [7[ but for globular clusters around 
NGC 3311 l lWehner et al j200a) . The model assumes a star forma- 
tion efficiency of /* = 0.2 and a pre-enrichment level of [m/H] = 
— 1.2 for the left sequence, with all other parameters identical to 
the model in Figure [7] 

or to mass loss experienced by these clusters since their 
formation, which would result in an overestimate of their 
mass today relative to their initial Mqc- 

An entirely new feature of our model is that it predicts 
the existence of a MMR for the metal-rich sequence of 
clusters. Because their pre-enrichment level is ~ 1 dex 
higher, the added factor of self-enrichment has a smaller 
amplitude than along the blue sequence and cuts in at a 
slightly higher luminosity, but in principle it should be 
observable. 

To date, we do not have much observational material 
for GCSs where the red sequence is seen to extend to high 
enough luminosities for its predicted MMR to show up. 
However, one such galaxy where hints of this effect can be 
seen is the Hydra cD NGC 3311, shown in Figure [5] The 
color-magnitude data (|Wehner et al. 1 120081 ) in (i 1 , g' — i') 
show the usual presence of the blue and red sequences, 
but the red sequence displays an unusually high upward 
extension that may co nnect to the UCD ma ss regime 
at 10 7 M© and above (jWehner fe Harris! [2007h and also 
shows a small 'swing' in mean color further toward the 
red. We have overplotted our model predictions, assum- 
ing the following transformations, 

%' = 37.90 -2.5 log— , (33) 

g'-i' = 1.173 + 0.311 log— (34) 

(based on lPeng et al.ll200"rl ). and assuming a star forma- 
tion efficiency /* = 0.2. The slight reduction in /* from 
the fiducial value of /* = 0.3 is required in order to fit 
the luminosity at which the MMR begins; the sensitive 
dependence of Mgp™ on /» in equation ([50]) means that 
very little change is required. In Figure [SJ it can be seen 
that the strong MMR of the metal-poor population and 
the mild MMR of the metal-rich population are both well 
reproduced in our model. We believe that new searches 
for this MMR effect at the high-mass end of both red and 
blue sequences would be an extremely effective test of the 
basic features of our model. They would best be looked 
for in individual supergiant ellipticals with the largest 
possible GC samples. 

6. CAVEATS 



6.1. Sensitivity to model parameters 

Our numerical results may be influenced by the par- 
ticular values we have adopted for parameters that are 
either unconstrained or that have significant uncertainty. 
We address here the influence that each of these param- 
eters has on our numerical and qualitative conclusions. 

6.1.1. Initial Mass Function 

The IMF is defined by three parameters: M m ; n , M max , 
and a. Changes of 20% in M min result in changes of 
5-10% in both the metallicity and its cluster-to-cluster 
scatter, while changes of 20% in M max result in changes 
of ~ 15% to the metallicity and its scatter. 

Flattening the slope of the IMF to a = —2.1 increases 
the mean metallicity by a factor of two, while steepening 
it to a — —2.5 decreases the mean metallicity by 35%; 
however, the changes to the scatter are in both cases less 
than 25%. Because changes to a result in changes to the 
mean metallicity that are comparable to the observed 
scatter, scatter in a between clusters could increase the 
predicted metallicity scatter to the observed level, as- 
suming that the other parameters of the IMF are kept 
constant. However, this is likely a poor assumption: the 
fraction of the stellar mass contained in high-mass stars, 
which determines the number of supernovae and there- 
fore the Poisson scatter, is a more robust quantity than 
the value of each parameter of the IMF. We conclude 
that our qualitative results are robust to any realistic 
uncertainty in the IMF. 

6.1.2. Supernova yields 

Because our equation for the hea vy-element yields is 
fit to a sm all number of models in IWooslev fe Weaver] 
(|1995() and iNomoto et all (|1997|), there is considerable 
uncertainty in the parameters. However, our results 
are completely insensitive to the B parameter in equa- 
tion ([3]), and while changes in the C parameter result in 
proportional changes to the mean predicted metallicity, 
they have no effect on the metallicity scatter. Another 
method of effectively changing the supernova yield is to 
change the minimum supernova progenitor mass, Msn- 
However, most of the metals are produced by stars with 
masses M* ^ Msn, and so changes in Msn have no not- 
icable effect on either the mean metallicity or its scatter. 

Another potential source of metallicity scatter is from 
metallicity-dependent supernova yields. If the mass of 
metals ejected from supernovae depended strongly on the 
metallicity of the progenitor star, then a small scatter 
in pre-enrichment levels could be leveraged into a large 
scatter in final metallicity. However, while the yields 
of individual elements vary significantly with progenitor 
metallicity, the total heavy elemen t yield s in the super- 
nova models of IWooslev fc Weaverl (|l995l) are very simi- 
lar for models with initial metallicities ranging from 10 -4 
to 10~ 2 Zq, ruling out this possibility. 

6.1.3. Star formation efficiency 

The predicted metallicity scatter for self enrichment is 
entirely independent of the mean value of the star for- 
mation efficiency, /*. Therefore, uncertainty in its value 
does not alter our conclusions. However, cloud-to-cloud 
scatter in the star formation e fficiency translates directl y 
into metallicity scatter (e.g. iRecchi fc Danzigerl [2005). 
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We therefore cannot rule out a large mass-independent 
variation in /* as an explanation for the observed metal- 
licity scatter among metal-poor Galactic GCs. However, 
our metal retention model of § [5] suggests that if self- 
enrichment were important, there would be an observ- 
able MMR for these clusters (only clusters with masses 
M > Mg£ aiu would exhibit no MMR, a regime that 
Galactic GCs clearly do not inhabit), which is not ob- 
served. 

6.1.4. Metal retention efficiency 

The metal retention efficiency, fz, occupies a similar 
role as /»: its mean value has no effect on the pre- 
dicted metallicity scatter in the self enrichment model, 
but cloud-to-cloud scatter in fz translates directly into 
metallicity scatter. We therefore cannot rule out a large 
mass-independent variation in fz as an explanation for 
the observed metallicity scatter, but again note that 
the MMR that would be expected in the case of self- 
enrichment is not observed for the Galactic GCs. 

Possible scatter in fz due to stochastic effects is dis- 
cussed in § I6.3( while scatter in fz due to scatter in 
cluster radii at a given mass is discussed in more detail 
in § 16.61 

6.2. Stellar winds 

We have assumed that the global heavy element en- 
richment and feedback energetics are dominated by core- 
collapse supernovae (i.e. Types lb and II). Stellar winds 
from stars of a wider range of mass contribute signifi- 
cant amounts of some individual elements, but are not 
believed to dominate the global enrichment and are cer- 
tainly unim portant energetically in comp arison to super- 
novae (e.g. lOppenheimer &: Davll2008h . If the metal 
contribution from stellar winds, i.e. from more com- 
mon lower-mass stars, is higher than currently believed, 
the metallicity scatter due to self enrichment would be 
further reduced from the values we derive. Our conclu- 
sion that self enrichment produces too little metallicity 
spread is therefore robust to the heavy element contri- 
bution from stellar winds. 

6.3. Stochastic metal retention 

In the previous picture, all supernovae simultaneously 
enrich and deposit energy into the protocluster cloud, 
which becomes well mixed and subsequently loses some 
fraction of its mass, and so a fraction fz of each super- 
nova's metals are retained. Another possible source of 
scatter is stochastic metal retention. In this picture, the 
individual supernovae pepper the cloud one at a time, 
and initially all metals are retained, but eventually one 
supernova injects enough energy that the total remain- 
ing gas can be unbound. Then, without any remaining 
gaseous envelope, the ejecta from all subsequent super- 
novae can escape freely and the enrichment stops. In 
this picture, a fraction fz of the supernovae have all of 
their metals retained, while no metals are retained from 
the remainder. The effective number of supernovae is 
then reduced from A^Nto /z-ZVsn, and the stochastic rel- 
ative scatter is increased by a factor of f z 1 ' 2 (this can 
be thought of as a special case of scatter in fz, which is 
discussed above). For fz — 0.08, the metallicity scatter 
would be increased by a factor of 3.5. This factor is well 



below the order of magnitude required to reconcile the 
observed and predicted metallicity scatter. 

6.4. Protocluster density profile 

The adopted isothermal density profile for the proto- 
cluster cloud, equation (|2Tj) , is as steep a profile as can 
be realistically maintained, and is physically well moti- 
vated. However, clouds may exist with slightly shallower 
density profiles. We therefore briefly expand the deriva- 
tion of § [5] to generic power-law profiles. 

If the density within the truncation radius is given by 

M c (3-/3) fry* . 
then for (3^2 the metal retention efficiency is 



fz 



1 



N SN E SN r t (2-(3) 
GM3 



(36) 



We have plotted this relation for a few slopes in Fig- 
ure 03 in addition to the fiducial (3 — 2 TSIS case. It 
is apparent that although the (3 = 2 case has a unique 
functional form, it is not "critical" in the sense of show- 
ing fundamentally different behaviour; in all cases, the 
metal retention turns over sharply near M Iotam . The 
qualitative effect of a shallower density profile is that a 
greater fraction of the gas is located near the edge of the 
cloud where it can be more easily unbound; fz there- 
fore decreases when (3 is lowered, and shows a sharper 
transition to the self-enriched regime. 

6.5. Mass-radius relation 

In § [5] we assumed that all protocluster clouds have 
a constant radius independent of mass. While this is 
true for observed GCs with M < 2 x 10 6 M Q , a mass- 
radius relation emerges at higher masses, with radii ris- 
ing from rh — 3 p c at ~ 10 6 Mp to r^ ~ 20 pc at 



, uum i h — o p c au ~ iu 1V1 G) m [ h — ^ pc au 
3 x 10 7 M M (see iReikuba et al.l l2007t iBarmbv et all 



l2007t lEvstigneeva et al.ll2008f l Assuming that the trun- 
cation radius scales similarly as the half-mass radius, we 
parametrize the mass-radius relation as follows: 

f r t , if Mgc < M lad , 

rt= l^fe)" ifM GC >M rad , ^ 

and examine the consequences on our previous results. 
The case n = 1/2 corresponds to simple virial equilib- 
rium and n = 1/3 to constant density. The available 
data, accumulated from a combination of the most mas- 
sive known GCs, UCDs, and dE nuclei, indicate n ~ 0.5 
(see the references cited above), which we adopt here. 
With this scaling, the metal retention efficiency be- 



fz = 



exp 



exp 



Air 



M GC 



if M gc < M rad , 



if M G c > M 



(38) 



rad i 



where we define the critical metal retention mass in the 
absence of a mass-radius relation M^f™ as in equa- 
tion (f30|) . The correct critical retention mass M retam 
is unchanged if Af ra( j > Mqqq*, but otherwise increases 
by a factor of 



M rotain 

l\ /"retain 
1V1 GC,0 



i\/f retain 
IV1 GC,0 



M r 



.1.1 



(39) 
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Fig. 9. — (a) Metal retention efficiency, fz, as a functio n of 
cluster mass, assuming the mass-radius relation of equation (1376 
(solid line). The dotted line indicates the relationship for clusters 
of constant radius, as in Figure [5] (b) Predicted mass-metallicity 
relation, as in Figure [6] but assuming the mass-radius relation of 
equation H37I I. The self enrichment level for clusters of constant 
radius is shown as the lower dashed line. 

For n = 1/2, the exponent in the equation above is just 
1. For example, if we adopt M rat j = 2 x 10 6 M as the 
transition to the mass-radius relation, then the critical 
retention mass increases by a factor of ~ 5 to Mg^, am = 
6 x 10 7 M Q . 

We have plotted our predicted fz and mass-metallicity 
relations in Figure [51 assuming the above mass-radius re- 
lation with n — 0.5 and M rad = 2 x 10 6 M Q . We have 
also plotted the original relations, where the cluster ra- 
dius is assumed to be constant, for comparison. We note 
that the increased radius of high-mass clusters reduces 
their metal retention because of the shallower potential 
well, and therefore moves the scale of self-enrichment to 
higher mass. However, it also results in a shallower mass 
dependence, more similar to the observed MMRs, and 
therefore the scale at which self-enrichment becomes ap- 
parent (which is approximately where the dotted and 
dashed lines cross in FigureEb) is less drastically affected 
than the factor of 5 difference in Mg^i ain . 

6.6. Scatter in radius 

Lastly, we consider the scatter in cloud radius at a 
given mass and its effect on fz- The observed scatter 
in r/j for GCs in the M ilky Way is ±0.3 in logr^ (with 
data from iHarrisl fl996i) : i.e., a r /r ~ 0.3. As discussed 
in § I6.1.4[ fz scatter translates directly into metallicity 
scatter. We therefore must determine whether the scat- 
ter in r t can help reconcile the self-enrichment model 
with the observed metallicity scatter. 

First, we note that the exponential dependence of fz 
on r t in equation (|28p implies: 



fz 



<7 r 

n 



ln/ z 



(40) 



From § 13.11 we determined fz < 0.08. With a r /r t ~ 
0.3, the implied scatter on fz is &f z / fz }t 0.75. This 
increase in predicted scatter is sufficient to account for 
the observed metallicity scatter amongst low mass GCs. 

While this appears to invalidate our argument in § 13.21 
that the magnitude of the observed metallicity scatter 
rules out self-enrichment as the dominant mode of metal 
enrichment, there are critical problems with this expla- 
nation. First, it is sensitive to the detailed form of fz, 
which we do not claim is robust. F or example, some 
previ ous studies of self-enrichment (e.g. lParmentier et al.1 
1999) have modelled the effects of supernovae as a step- 
function switch where the protoclustcr is completely dis- 
rupted below a critical mass, but survives intact with 
full metal retention (fz = 1) above it; in such a case, 
the scatter in r t would not increase the metallicity scat- 
ter of the extant clusters at all. Second, the increased 



scatter is a direct result of the fact that fz is a function 
of the depth of the potential well, which is proportional 
to M c /r t . In other words, it is intimately tied to the 
dependence of fz on M c and the existence of a mass- 
metallicity relation; however, no mass-metallicity rela- 
tion is observed among the low-mass globular clusters. 
The only possible ways that this could be masked is if 
M c oc r t (in marked contrast to the observations, which 
indicate no mass-radius relation in this mass regime) , or 
if there were a fortuitous cancellation between /„ and fz 
of the form oc f% 1 , which we consider highly unlikely 
We conclude that the overall picture remains most 
consistent with metallicities that are dominated by pre- 
enrichment for most GCs. 

7. DISCUSSION 

7.1. Comparisons with Previous Models 

Several previous models for self enrichment have been 
propo s ed in the litera t ure (lCavrellll986l: [Morgan fc Lake! 



1989b iBrowneta l 119911 : IParmentier et I all Il999t 



Recchi fc Danzigerll2005UStrader fc Smith! l2008f ) Most 
of these are based on the picture of iFall fc Reesl (|1985h . 
in which the protocluster cloud is a pressure-truncated 
condensation in the hot halo rath er than having its 
structure determined by gravity (iBrown et al. 1991; 
IParmentier et al.l Il999b IRecchi fc Danzigeri I2005T) . Such 
models have pred icted no mass-metallicity relation 
(|Brown et al.l ll991). or a minimal blue tilt that is m uch 
weaker than that observed (|Recchi fc Dan zigcr 20QH), or 
an inverse MMR in which lower-m ass clusters are more 
metal- rich (jParmentier et al1ll999h . These are therefore 
not adequate for understanding the metallicities of the 
high-mass clusters whose properties are explained by our 
model. Most of these models also invoke two distinct 
star formation episodes: an initial burst of (possibly 
massive) primordial composition stars whose ejecta 
pollute the protocluster region, followed by a subsequent 
burst of star formation from the enriched material. We 
would classify such a model, in which the metals are 
produced prior to the burst that creates the stellar 
cluster we see today, as "pre-enrichment" ; however, if 
the delay between the two bursts is sufficiently short 
(of order the dynamical time), then this is a m atter of 
semantics rather than a physi cal difference (e.g. ICavrell 
[19861: IRecchi fc Danziged^OOl. 

In the model of IRecchi fc Danzigerl (|2005l ). they sug- 
gest a number of alternatives in order to obtain the right 
cluster-to-cluster metallicity scatter er[Fe/H], the most 
notable of which is to adopt arbitrary and large cluster- 
to-cluster differences in their "mixing efficiency" , a pa- 
rameter resembling our gas retention fraction fz or the 
effective yield y e f f. If fz is primarily determined by the 
initial cluster mass M c , as we suggest that it should, then 
we view this approach as unlikely. Instead, if er[Fe/H] is 
due to pre-enrichment, then it can be understood classi- 
cally as the result of simple chemical evolution such as 
a "leaky box" model with a suitable effective yield. For 
the metal-poor GCs or the Milky Way halo stars, a rea- 



sonable choice is y e jf ~ 0.03 (Rvan fc Norrislll991l : 



Prantzos 2003; Vandalsen & Harris 



2004). 



For o ur purposes, the most re levant previous model is 
that of IStrader fc S mith (2008), which was specifically 
created to explain the observed MMR. Their model, like 
ours, invokes energetic balance between supernova feed- 
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back and gravity within a single episode of star forma- 
tion to produce an MMR via self enrichment. They also 
put forward the ideas that the transition mass marks the 
point at which the metals produced by self-enrichment 
become comparable to those due to pre-enrichment, and 
that the lack of a corresponding "red tilt" is due to the 
higher level of pre-enrichment among red sequence GCs. 
However, there are several key differences between our 
models. First, in their model the star formation efficiency 
/* is assumed to vary with mass while the metal retention 
fz is constant; instead, we argue that fz is the physical 
quantity much more likely to depend strongly on clus- 
ter ma ss, whereas /* ~ const. Second, IStrader fc Smithl 
(2008]) derive power-law scalings for the cluster proper- 
ties such as mass, metallicity, and radius, but do not ad- 
dress their amplitude at all; in contrast, our model is on 
an absolute scale and explicitly predicts the transition 
mass from the pre-enriched regime to the self-enriched 
regime. Third, we predict the existence of a more mod- 
est MMR at the high-mass end of the red GC sequence. 
Finally, and most importantly, they examine the mean 
trends but not the metallicity scatter between clusters; 
not only do we explicitly predict the metallicity scat- 
ter due to self enrichment, but we use this prediction 
to conclude that self-enrichment is unimportant among 
low-mass GCs. 

It may be that a combination of ideas from different 
models would produce an even more realistic match to 
the real GC distributions. For example, our model pre- 
dicts a relatively rapid transition from the pre-enriched 
regime at low cluster masses to a strongly self-enriched 
regime at high masses, which may be more rapid than 
observed (note, however, that the mass-radius relation 
described in § 16.51 results in a shallower MMR) . How- 
ever, a rapid transition from no metal retention to sig- 
nificant metal retention, as in our model, coupled with 
a star formation efficie ncy that scales with mass, as in 
IStrader fc Smithl J2008), might result in a model that is 
in even better agreement with the observations. 

Another possible feature to include in a more general 
model might be the role of external pressure confinement, 
which could help the lower-mass clouds keep back some of 
their SN ejecta (rather than having fz — ► as would be 
the case if only self-gravity were operating) and possibly 
allowing the MMR to extend to somewhat lower mass. 

7.2. Variation Between Galaxies 

At present, the different sets of observations indi- 
cate that the amplitude (i.e. the exponent in the 
scaling Zqc °c Mgc") of the MMR may differ from 
galaxy to galaxy. The possibility also exists that 
the blue-sequence MMR may set in at different mass 
scales and in som e cases may be completely absent 
(lHarr s et alJ 120061: iMieske et al.l 120061: IStrader et all 



120061: iSpitler et all [20061 120081: iWehner et all l2008TT ~Ln 
the context of our simple enrichment-based model, these 
differences would boil down to differences in the degree 
of self-enrichment at a given cluster mass. The following 
factors could vary from galaxy to galaxy and therefore 
affect the observed cluster metallicities: 



present-day masses less resembling their original 
Mqc- In this case, the MMR imparted at birth 
would be weakened or destroyed when compar- 
ing the present-day masses of the clusters to their 
metallicity. A possible consequence of this effect 
would be that the MMR should be stronger in the 
outer halos of large galaxies. 

2. If the protoclusters were bathed in a particularly 
strong ultraviolet background during formation, 
(e.g. if the galaxy was undergoing a massive star- 
burst during the main epoch of cluster formation) , 
the gas may have been less able to cool and there- 
fore the star formation efficiency /* may have been 
reduced. This would result in a lower critical mass 
M retam and a lower Z c . If /* were reduced suffi- 
ciently, the maximum metallicity achievable by self 
enrichment could even drop below the pre-enriched 
value, resulting in no MMR at all (see also the dis- 
cussion about the influ ence of a UV background in 
IStrader fc Srmthll2008h . 

8. CONCLUSIONS 

We have developed a model for stochastic self enrich- 
ment in globular clusters. This model predicts both 
the mean metallicity and the cluster-to-cluster spread in 
metallicities that is expected due to stochastic sampling 
of the IMF. The predicted metallicity scatter is an order 
of magnitude smaller than observed for Galactic GCs; 
this rules out self enrichment as an important contributor 
to their global metal content and leaves pre-enrichment 
as the dominant contributor. This conclusion does not 
depend strongly on the adopted value of any free param- 
eter and is robust to reasonable changes to the shape of 
the IMF, the supernova heavy element yields, and the 
details of how metals are retained. Although significant 
cluster-to-cluster scatter in the star formation efficiency 
and/or the fraction of metals retained at a given cluster 
mass could increase the predicted scatter to the observed 
level, they would most likely be accompanied by a sig- 
nificant mass-metallicity relation at relatively low mass, 
which is not observed for the Milky Way. 

We have used simple energetics to predict how the 
metal retention efficiency, a key parameter in the self- 
enrichment model, increases with cluster mass. We 
propose a combined pre-enrichment plus self-enrichment 
model, where clusters are pre-enriched to the observed 
level with significant scatter due to their environment, 
and then those clusters sufficiently massive that they can 
retain a significant amount of their supernova ejecta are 
further self-enriched; the threshold for measurable self- 
enrichment is predicted to be typically a few 10 6 Mq. 
This model matches the main features of the data so far. 

The key features of our model can be summarized as 
follows: 

1. For globular cluster masses less than ~ 10 6 Mq, no 
MMR is expected for either the blue or red clusters; 
the mean cluster metallicities should not depend on 
mass. 



1. Mass loss from a cluster depends on the strength 
of the tidal field in which it orbits. Therefore, 
GCs subjected to stronger tidal fields should have 



2. Along the blue GC sequence, the MMR should 
'kick in' noticeably for Mqc higher than a few 
million Mq, reaching [m/H] ~ —1 at A/qc ~ 
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1 — 2 x 10 7 Mq, corresponding to the most massive 
known GCs. 

3. The cluster-to-cluster metallicity spread tr[m/H] 
remains uniform in the low-mass, pre-enrichmcnt 
regime, but decreases at higher mass in the MMR 
regime where self-enrichment becomes most impor- 
tant. 

4. A more modest MMR should exist for the red GC 
sequence, beginning at a few million Mq. In at 
least one cD galaxy with a very rich GC system 
(NGC 3311), this effect at the top of the red se- 
quence may have already been seen. 

5. Within the context of our model, the main free pa- 
rameter controlling the level and amplitude of the 
MMR is the star formation efficiency /*. Lower- 
ing /* from our baseline value of 0.3 would reduce 



the slope of the MMR along both the red and blue 
GC sequences and might explain the observed dif- 
ferences between galaxies. 

Two other physical effects not included in our basic 
model are a possible correlation of star formation effi- 
ciency with cluster mass, and pressure confinement from 
the gas outside the protoclusters at time of formation. 
Including these effects properly would allow for a greater 
range of MMR slopes and onset points and would be an 
obvious route to explore in a next stage of studying this 
intriguing effect. 
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